Fullerene-based molecular nanobridges: A first-principles study 
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Building upon traditional quantum chemistry calculations, we have implemented an ab-initio 
method to study the electrical transport in nanocontacts. We illustrate our technique calculat- 
ing the conductance of Ceo molecules connected in various ways to Al electrodes characterized at 
the atomic level. Central to a correct estimate of the electrical current is a precise knowledge of the 
local charge transfer between molecule and metal which, in turn, guarantees the correct positioning 
of the Fermi level with respect to the molecular orbitals. Contrary to our expectations, ballistic 
transport seems to occur in this system. 



At the forefront of the carbon-based molecules that 
are expected to play a key role in nanoelectronic devices 
one can situate carbon nanotubes and fullerenes. From 
Coulomb blockade to Kondo-like phenomena through 
ballistic transport, multiple transport regimes have been 
observed in nanotubes (even for the same chirality of the 
nanotube) depending on the way the contact with the 
electrodes is made Q. Very little is known, however, 
about the true role played by this contact as well as of 
the effect of the substrate on which the carbon nanotube 
has been layed. This is somewhat connected with the 
difficulties experienced in the interpretation of electronic 
transport experiments through Ceo molecules. Experi- 
ments using a scanning tunneling microscope (STM) Q 
have shown the possibility of "looking" at individual Ceo 
molecules when they are adsorbed on a conducting sub- 
strate. There is still, however, a lot of controversy in the 
interpretation of the images obtained by different groups 
H. Alternative to STM set-ups, nanoscopic break junc- 
tions also revealed themselves as powerful tools to study 
transport through individual molecules f|,||. Recently, 
for instance, in between the two electrodes of a gold 
break junction, evaporated Ceo molecules unexpectedly 
created "mechanical" bridges for electrical transport be- 
tween electrodes |J. 

Archetypal STM and break-junction set-ups above 
mentioned can be cataloged under the term molecular 
nanobridge, i.e., metallic electrode + molecule + metallic 
electrode. The conductance, G, is essentially determined 
by the electronic structure of the isolated molecule since 
the number of available channels for conduction deep in 
the electrodes is always much larger than those provided 
by the molecule. However, as already mentioned, this 
conductance is strongly influenced by the chemistry of 
the electrode-molecule contacts. Behind this contact lies 
a key problem: Charge transfers between electrode and 
molecule when they come in close proximity. The ne- 
cessity of aligning the Fermi level of the metallic elec- 



trodes with the "Fermi level" of the molecule forces a 
certain amount of charge to be transferred one way or 
the other. The way this alignment is achieved and how 
much charge is transferred constitutes a difficult prob- 
lem which turns out to be essential to understand the 
electrical transport in molecular nanobridges and, more 
generally, in nanocontacts or nanoconstrictions. Curi- 
ously enough, the prevailing theoretical descriptions in 
the great body of work done on various types of metal- 
molecule-metal systems |?J as well as on atomic chains 
and constrictions || are still based on parametrized tight- 
binding or semi-empirical models. The theoretical frame- 
work to calculate the current is well established [[)), but, 
within these models, Fermi level alignment has to be ei- 
ther entirely ignored or imposed with some additional 
criteria. While local charge neutrality is a sensible cri- 
terion for atomic chains and constrictions made of the 
same material as the electrodes ||, there is not the like 
for metal-molecule-metal systems. This is, in part, the 
motivation for developing ab-initio methods which can 
provide qualitative and quantitative answers in all pos- 
sible scenarios. Some groups have already blazed the 
trail in this direction pp|-|l^|, although in their methods 
the problem of the atomic structure of the electrodes is 
usually put aside. This is not satisfactory when one is 
trying to describe STM experiments where the detailed 
atomic structure of the tip determines, to a large extent, 
whether or not the STM images can resolve the topogra- 
phy or molecular structure of the adsorbate. 

What we present here is an ab-initio alternative to 
parametrized tight-binding or semi-empirical models in 
order to address electrical transport through generic 
nanocontacts. Our method automatically accounts for 
Fermi level alignment and charge transfer whenever the 
latter is expected to occur. This is the case in the 
system we study here: Aloo-Ceo-Aloo- Isolated Ceo 
molecules present a large gap (~ 2 eV.) and behave like 
intrinsic semiconductors. At low temperatures, metal- 
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semiconductor-metal systems do not conduct when the 
semiconductor is intrinsic regardless of the sign and the 
amount of charge transferred at the interfaces. In the 
met al-fullerene- metal system, however, a tunnel current 
can always be expected due to the nanoscopic size of 
the fullerene JH| ]. What is more surprising is the fact 
that Ceo molecules can conduct ballistically (G > e 2 /h) 
. Some of the lowest unoccupied molecular orbitals 
(LUMO's) of the isolated fullerene still retain its ex- 
tended character when contact is made to the electrodes 
and they become partially or fully occupied. In solid- 
state language: The Fermi level lies above the bottom 
of the conduction band all across the (nanoscopic) semi- 
conductor. Notice that bigger molecules exhibiting a gap 
like some families of undoped nanotubes are not expected 
to conduct at zero bias [Q. Whether or not this occurs 
depends on the size of the molecule and the details of 
the charge transfer process which, in turn, relies heavily 
on the relaxation of the atomic structure of the fullerene 
when in contact with the metal. 




FIG. 1. Schematic view of the two-step process for calculat- 
ing the Green's function used in the calculation of the conduc- 
tance: (a) A density functional calculation of a Ceo molecule 
contacted by two electrodes (here in the form of pyramids) is 
performed, (b) Atoms that are irrelevant for electrical trans- 
port are removed and substituted by an effective selfenergy 
"attached" to the contact left and right atoms remaining in 
the electrodes (see text). 

Standard ab-initio or quantum chemistry calculations 
[|l5| , p!6[ are only possible for finite or periodic systems, 
whereas the transport problem one typically wants to ad- 
dress requires infinitely large leads with no symmetry at 
all. Drawing on Refs. [|l2],[l3|, the main idea here consists 
of performing a density functional (DF) calculation of the 
molecule including part of the leads with the desired ge- 
ometry see Fig. |(a)] . This can be efficiently done using 
the Gaussian-98 code Gq|. For a sufficiently large num- 



ber of atoms describing the metallic electrodes the en- 
ergy of the highest occupied molecular orbital (HOMO) 
of the whole structure sets the Fermi level. In addition, 
the correct charge transfer (to or from the molecule) is 
guaranteed. Next, we perform a Lowdin orthogonaliza- 
tion of the original non-orthogonal gaussian basis set into 
an orthogonal one which, having basis elements more ex- 
tended than the original one, still preserves the symmetry 
of the orbitals and the atomic character. The hamilto- 
nian of the electrode+molecule+electrode system, H, as 
it stands, is finite and, according to the usual theoreti- 
cal transport schemes ||, its Green's functions are un- 
suitable for any current determination since they simply 
have poles. The following (and crucial) step in our proce- 
dure is to transform this finite system into an effectively 
infinite one. In order to do this we first remove from the 
hamiltonian all but the N atoms forming the relevant 
atomic structure of each electrode close to the molecule 
see Fig. 0(b)]. The retarded Green's function associated 
with this reduced hamiltonian H r is now transformed 
into the retarded Green's function of an infinite system: 

G r (e) = (el - H r + id)' 1 -> [el - H r - E(e)]- 1 . (1) 

In this expression £ = E# + E^ where Y>r(Y>l) denotes 
a self-energy matrix that accounts for the right (left) in- 
finite electrode, part of which has been included in the 
DF calculation [see Fig. 0(a)]. The added self-energy can 
only be explicitly calculated in ideal situations, which, in 
principle, might seem to limit the desired applicability 
of the above step. However, as pointed out by Landauer 
in his scattering description of electrical transport [[)), 
the bulk details of a metallic electrode are not neces- 
sarily relevant for the electrical transport properties of 
the region with the smallest number of channels, i.e., 
the contacted molecule. With this crucial observation 
in mind we choose to describe the bulk electrode with 
a Bethe lattice tight-binding model JTt| with the appro- 
priate coordination and appropriate parameters. More 
specifically: We require our Bethe lattice model to re- 
produce the electrode bulk density of states and to have 
the same Fermi energy as that of the system on which 
the DF calculation was initially performed |l8| . The ad- 
vantage of choosing a Bethe lattice resides in that E can 
be easily calculated through a well-known iteration pro- 
cedure JTt| which we do not detail here. It suffices to 
say that, for each atom in the reduced electrode that has 
been stripped from, at least, one nearest-neighbor atom, 
a contribution from an infinite Bethe lattice calculated in 
the direction of the removed atom is added to E. Assum- 
ing that the initial system was large enough and that the 
reduced system still contains the relevant atomic details 
of the electrode close to the contact with the molecule, 
to replace a piece of the electrode by our effective self- 
energy should introduce no spurious effects, and it makes 
the system effectively infinite. The conductance can now 
be simply calculated through the expression B 
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where Tr denotes the trace over all the orbitals of the 
reduced hamiltonian and the matrices Tr and Tl are 
given by z(Xl^ — S^) and i(Tt r L — respectively. As 
in the self-energy matrices, the matrix elements of Tr 
and are only different from zero for the atoms in the 
reduced electrodes that have been stripped from at least 
one near-neighbor atom. 
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FIG. 2. (a) Conductance around the Fermi energy (here set 
to zero) of a Ceo molecule contacted by two opposite Al elec- 
trodes in the form of pyramids composed of 30 atoms each 
(see inset). We present three different cases: N = 1 (dot- 
ted line), N = 5 (dashed line), and N = 14 contact atoms 
(solid line), corresponding to an increasingly detailed elec- 
trode structure in the reduced system (see text), (b) Same as 
in (a), but with the fullerene in between Al(001) parallel sur- 
faces, each one composed by two atomic planes with a total 
of 33 atoms (see inset). We present results for N = 1 (dotted 
line), N = 5 (dashed line), and N = 9 contact atoms (solid 
line), corresponding to an increasing surface plane. 

We apply here the technique just described to the sys- 
tem Aloo-Ceo-Aloo with different geometries. Although 
they are not commonly used in experiments, we have 
considered Al electrodes for two independent reasons: (i) 
The number of electrons per atom is small enough to 
avoid pseudo-potentials in the DF calculations and, (ii) in 
contrast to noble metals, Al presents a strong, essentially 



covalent chemical bond with carbon molecules which, in 
principle, makes it ideal for the future fabrication of sta- 
ble nanostructures. For the DF calculation we have used 
the Becke's three-parameter hybrid functional using the 
Lee, Yang and Parr correlation functional with a mini- 
mal basis set [|l9[ . The orientation of the fullerene and the 
separation and atomic structure of the electrodes consid- 
ered in Fig. ^ (see insets) are reasonable, but arbitrary 
[ pQ[ . In Fig. |(a) we plot the conductance around the 
Fermi energy (here set to zero) of the fullerene contacted 
by two opposite electrodes in the form of pyramids [see 
inset or Fig.|l|(a)]. These are composed of 30 atoms each 
and their apex atoms are separated by 10.1 A . We have 
orientated the fullerene so that two opposite C atoms 
are as close as possible to the apex atoms of the pyra- 
mids (i.e., 1.4 A between the C atom of the fullerene and 
the Al apex atom). In order to illustrate the reliability 
of the method, we present the results for three differ- 
ent approximation degrees: N = 1 (dotted line), N = 5 
(dashed line), and N = 14 (solid line), corresponding to 
an increasingly detailed contact structure. Notice how, 
as expected, all the curves are very similar since the rel- 
evant atomic structure close to the molecule is simply 
represented by the atom at the pyramid apex. Logically, 
the results using 5 atoms (two atomic planes) are almost 
identical to those using 14 atoms (three atomic planes). 
The gap of the isolated fullerene reflects itself in a van- 
ishing G below the Fermi energy and the LUMO's of the 
isolated molecule are now partially occupied. Note that 
they retain their extended character since the conduc- 
tance around the Fermi level goes up to 2e 2 /h at cer- 
tain energies. Nonetheless, the three-fold degeneracy of 
the fi u and f\ Q LUMO's has been entirely removed by 
the interaction with the electrodes (this interaction is 
also responsible for the broadening of the peaks). In this 
situation Coulomb blockade physics must dominate and 
spin-splitting will occur in the fullerene [^l| . 

In Fig. ||(b) we present the conductance of the fullerene 
oriented as in the previous case, but where the two pyra- 
mids have been replaced by Al(001) parallel surfaces (see 
inset). The different types of lines correspond to an in- 
creasing definition of the surface plane (1,5, and 9 surface 
atoms). Now one atom is not enough to define the con- 
tact since there is an important coupling between several 
Al surface atoms and the fullerene. Not even 5 atoms suf- 
fice to define the metallic contact region quite precisely. 
Now the value of the conductance rises up to « 2.5e 2 /h 
right at the Fermi energy. This was somehow expected 
since, in this situation, the hybridization with the surface 
is stronger and there are more channels for the electron to 
travel between electrodes while the LUMO's appear to be 
also extended. The gap still reflects itself in G, but has 
been partially closed by the strong hybridization. The 
excess charge in the fullerene is ~ 2.6 which agrees fairly 
well with a related and recent calculation [01. 
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FIG. 3. Conductance around the Fermi energy (here set 
to zero) of a Ceo molecule contacted by two parallel Al(lll) 
surfaces separated by 10.1 A and composed of 31 atoms each 
(see inset). The dashed line is the result for the non-relaxed 
fullerene and the solid one is the result when the fullerene has 
been allowed to relax. 

Figure || shows the conductance when the fullerene is 
placed in between two parallel Al(lll) surfaces (31 atoms 
each, 7 of which have been used in the calculation of G) 
with parallel Al surface and fullerene hexagons. This 
orientation has been chosen on the basis of several ex- 
perimental observations In one case (dashed line) 
the fullerene has not been allowed to relax. In the other 
(solid line) an ab-initio relaxation has been performed. 
The magnitude of G in the former case is comparable 
to that shown in Fig. §(b). Again, a large number of 
surface and fullerene atoms intervene. While the de- 
tails change, the gap is still visible around -2 eV and 
the excess charge is « 2.1. Surprisingly enough, when 
the fullerene is allowed to relax (stretch in this case), 
the conductance changes drastically. Only one LUMO 
around the Fermi level preserves its extended character, 
the rest becoming localized at the interfaces. Both a Mul- 
liken population analysis and an integration of the local 
density of states show that the excess transferred charge 
(~ 3.3) concentrates now more at the interface hexagons 
where bonding takes place. (The distance between sur- 
faces was chosen to be slightly larger than the diameter 
of the fullerene so the structural changes of the fullerene 
do not play any role). The small value of the conduc- 
tance below the Fermi energy cannot be attributed to 
the gap of the fullerene (which has entirely disappeared 
from the density of states), but to the localization pro- 
cess which does not reflect in the global density of states 
of the fullerene. This last result bridges the gap with the 
intuitive picture for transport in a metal-semiconductor- 
metal system and shows that the chemical details of the 
metal-molecule contact are determinant in the transport 
properties of molecular bridges. 

In summary, we have developed a methodology to cal- 
culate ab-initio transport properties of nanoconstrictions. 



We have chosen to study a system of recent interest: A 
fullerene contacted by Al electrodes. Our main goal was 
to describe correctly the contact charge transfer so that 
realistic predictions on the electrical conduction proper- 
ties of these molecules can be made. 
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